Density Functional Study of the Over-Doped Iron Chalcogenide: TlFe 2 Se 2 with 

ThCr 2 Si 2 structure 
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We report density functional calculations of electronic structure and magnetic properties of ternary 
iron chalcogenide TlFe2Se2, which occurs in the ThCr2Si2 structure and discuss the results in rela- 
tion to the iron-based superconductors. The ground state is antiferromagnetic with checkerboard 
order and Fe moment ~ 1.90 [iB. There is strong magnetoelastic coupling similar to the Fe-based 
superconductors, reflected in a sensitivity of the Se position to magnetism. Tl is monovalent in this 
compound, providing heavy electron-doping of 0.5 additional carriers per Fe relative to the parent 
compounds of the Fe-based superconductors. Other than the change in electron count, the electronic 
structure is rather similar to those materials. In particular, the Fermi surface is closely related to 
those of the Fe-based superconductors, except that the electron cylinders are larger, and the hole 
sections are suppressed. This removes the tendency towards a spin density wave. 

PACS numbers: 74.25. Jb,71.20.Lp,74.20.Mn,75.10.Lp 



I. INTRODUCTION 

The finding of superconductivity in layered Fe oxy- 
pnictides with the ZrCuSiAs structure^ 2 - has stimulated 
much interest leading to the discovery of many addi- 
tional compounds with critical temperatures up to ~ 
56 Kj 3 -^ Superconductivity occurs quite generally in 
these compounds, and is robust against off-site sub- 
stitutions, and variations in the structure, including 
doped fluoro-arsenides SrFeAsF,— £ as well as ThCr2Si2- 
type arsenides (prototype BaFe2As2 7 ), Cu2Sb-structure 
LiFeAs 8 . (NaFeAsS), and PbO structure Fe(Se,Te),±&±i 

All these materials have square planar Fe 2+ layers, 
tetrahedrally coordinated by anions and are near mag- 
netism. The undoped compounds generally show a spin 
density wave (SDW) ) 12 : 13 ' 14 and superconductivity ap- 
pears when this SDW is suppressed either by pressure 
or doping. This suggests an involvement of magnetism 
in the pairing. 15 In addition there is accumulating ev- 
idence for strong antiferromagnetic correlations in this 
family in the normal metallic states, e.g. from the tem- 
perature dependence of the susceptibility, which increases 
with T^L 

As mentioned, the binary chalcogenide a-FeSe is also 
superconducting when doped by off-stoichiometry^ T c 
for this material reaches 27 K under pressure,— while the 
electronic structure is very similar to the pnictide phases, 
according to both density functional calculations^ 8 - and 
photoemission^ Among the chalcogenides, experiments 
indicated that FeTe might have the strongest supercon- 
ductivity based on the increase in T c upon alloying FeSe 
with Te, (increasing from 8 K to 15 K)2S as well as theo- 
retical considerations based on the strength of its SDW. 18 
However, FeTe shows an antiferromagnetic ground state 
instead of superconductivity ) 21 : 22 ' 23 which implies addi- 
tional doping (or pressure) is still required to further sup- 
press magnetism and perhaps induce superconductivity. 
Modest superconductivity (10 K) has already been ob- 
served when alloying FeTe with Si2£ In this regard, it is 



important to note that a-FeSe and FeTe generally form 
off stoichiometry with excess Fe partially filling intersti- 
tial sites in the chalcogen layer. This excess Fe acts as 
an electron dopant, and suppresses SDW order by do- 
nating carriers to Fe-Se layers^ In addition, the excess 
Fe, carries local magnetic moments that may be expected 
to interact with spin fluctuations associated with the Fe- 
Se layers and also produce pair breaking ] 22 : 25 : 26 There- 
fore, it is of interest to explore alternate chemical doping 
strategies for FeSe and FeTe, perhaps based on related 
compounds. 

Here we report density functional studies of TlFe2Se2, 
which is a ternary iron chalcogenide occurring in 
ThCr2Si2 structure. Little is known about the physical 
properties of this compound, except that it has been syn- 
thesized and that it has a magnetic transition at ^450 K, 
as observed in Mossbauer measurements. The magnetic 
state was characterized as antiferromagnetic based on the 
lack of attraction of the magnetic phase by a magnet *2Z 



II. STRUCTURE AND METHODS 

From a structural point of view, TlFe2Se2 (Fig. [1} is 
closely related to FeSe, and in particular consists of prac- 
tically identical Fe-Se layers, consisting of edge-sharing 
FeSe4 tetrahedra, although these are now intercalated 
with Tl, altering the stacking sequence. The large size of 
Tl cations (even larger than Ba in BaFe2As2) results in a 
larger separation (7.0 A) between Fe-Se layers than FeSe 
(5.5 A). In our calculations, we used the experimental 
values of tetragonal lattice constants, a = 3.89 A and c 
= 14.00 A, 28 while the internal coordinate, zs c was re- 
laxed by energy minimization (the atomic coordinates in 
this structure are Fe (4d) (0, 0.5, 0.25), Tl (2a) (0, 0, 0), 
and Se (4e) (0, 0, z Sc ) ) 

The present calculations were performed within the 
generalized gradient approximation (GGA) of Perdew, 
Burke, and Ernzerhof (PBE)^ using both the general 
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FIG. 1: (color online) Crystal structure of ThCr2Si2-type 
TlFe 2 Se 2 . 



potential linearized augmented plane-wave (LAPW) 30 i 31 
and projector augmented-wave (PAW) methods. 32 ' 33 For 
the LAPW method, we employed LAPW spheres of ra- 
dius 2.5ao for Tl and 2.1ao for Fe and Se. Well converged 
basis sets with the size determined by i?Fefcmax=8.0 were 
used and semicore states of 5d for Tl, 3p for Fe, and 3d, 
4s for Se were included using local orbitals. In the PAW 
calculations, a kinetic energy cutoff of 350 eV and aug- 
mentation charge cutoff of 511 eV were used. A 16x16x16 
grid was used for the body-centered tetragonal Brillouin 
zone sampling in the self-consistent calculations, while 
more dense grids were used for density of states (DOS), 
Fermi surface, and especially magnetic calculations. The 
electric-field gradient (EFG) calculations were done with 
the all-electron LAPW method. The formation energy 
for vacancies was calculated with the PAW method us- 
ing a 90 atom supercell, not including magnetism. The 
chemical potentials of bulk elemental Tl and Fe was used 
as a reference. For bec Fe our calculated PAW spin 
magnetic moment of 2.27 fiB is in reasonable agreement 
with experimental value 2.12 fiB. Consistency between 
the LAPW and PAW calculations was carefully cross- 
checked and is indicated for example by very small resid- 
ual forces for the relaxed structure, as well as very closely 
coincident band structures and density of states with two 
methods. 
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FIG. 2: (color online) Calculated electronic band structure of 
TlFe2Se2 with the checkerboard antiferromagnetic order, us- 
ing the spin-polarized relaxed internal coordinate z$ c ~ 0.348. 



III. GROUND STATE: CHECKERBOARD 
ANTIFERROMAGNETISM 

We begin with the magnetic order of stoichiomet- 
ric TlFe2Se2, which is after all a known magnetic ma- 
terial. The energetics were calculated with optimiza- 
tion of the internal coordinate zs c considering three 
different long-range magnetic orders possibly existing 
within Fe sheets; 15 ' 34 ferromagnetism, nearest-neighbor 
checkerboard antiferromagnetism, and stripe antiferro- 
magnetism (the pattern of the SDW in the undoped Fe 
superconductors). We found that the energy is always 
lowered in magnetic configurations where the Fe layer 
stacking along the c axis in antiferromagnetic, i.e. the in- 
terlayer Fe-Fe interaction is antiferromagnetic. Relative 
to the non-spin-polarized case, we found a ferromagnetic 
instability (-44 meV/Fe) at the relaxed zs c — 0.360, with 
Fe moment 2.79 /iB/Fc. However a stronger instability 
is found for the checkerboard nearest neighbor antifer- 
romagnetic order (-78 meV/Fe, 1.90 /iB/Fe) at relaxed 
^Sc = 0.348, which is the ground state according to our 
calculations. This seems to be consistent with the only 
available experimental result, i.e. antiferromagnetism. 27 
Obviously, further experiments will be helpful in confirm- 
ing the magnetic order in TlFe2Se2. Interestingly, we did 
not find a magnetic instability at all for the SDW-like or- 
dering. The sensitivity of the moment formation to the 
ordering pattern shows that the magnetism is itinerant 
in nature. In non-spin-polarized calculations, we obtain 
a relaxed optimal internal Se coordinate of zsc=0.342, 
lower than both the the experimental value 0.357 (with 
the difference of 0.2 A) and the value obtained in cal- 
culations including magnetism, although with the lowest 
energy checkerboard order, we still obtain a value sig- 
nificantly lower than experiment. In any case, consid- 
ering also the difference between the ferromagnetic and 
checkerboard values of zs c we do find the substantial 
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magnetoelastic coupling, similar to the Fe-based super- 
conducting phases. 

The calculated band structure and electronic density 
of states (DOS) with the ground state checkerboard an- 
tiferromagnetic order are shown in Figs. [2] and [3l re- 
spectively. This material is metallic with two dispersive 
electron bands crossing the Fermi level {Ep) around the 
X point. As may be seen, the Fermi energy occurs in 
a trough in the DOS with a sharp peak ~ 0.5 eV below 
Ep. This peak arises from the very flat band below Ep, 
as seen in Fig. [2l At the Fermi energy, N(Ef)=1.1 eV -1 
per Fe. The Fe 3d spin-down states are almost filled. 
Integration of partial spin states up to E-p and normal- 
ization with total Fe 3d states give 4.3 electrons in the 
majority states and 2.4 electrons in the minority states. 
This indicates the valence state of Fe is lower than +2 
here. The Se p states are mainly found below -3 eV rel- 
ative to the Ep, and are only modestly hybridized with 
Fe d states. The Tl 6s states are occupied, while the 
remaining Tl 6p states are above Ep, indicating that Tl 
occurs as monovalent Tl + . This is consistent with elec- 
tron doping of the the Fe-Se sheets, by 0.5 e / Fe, i.e. 
nominal Fe valence of Fe 15+ . 

To connect with future experiments, we calculated the 
electric field gradient (EFG) of Fe and Se sites. These 
arc defined as the second derivative of the Coulomb po- 
tential at the nuclear sites, and are probed by nuclear 
magnetic resonance (NMR) or nuclear quadrupole reso- 
nance (NQR) measurements. For tetragonal site sym- 
metry, only the V zz component is independent since the 
EFG tensor must be traceless. The calculated V zz with 
the ground state antiferromagnetic order are -2.2xl0 21 
V/m 2 and -5.6xl0 21 V/m 2 for Fe and Se, respectively. 
As noted below these values are very sensitive to the 
value of the internal coordinate, and therefore EFG mea- 
surements may be a very useful probe of magnetoelastic 
coupling in these materials. In this regard, it should 
be noted that Mukuda and co-workers observed a rela- 
tionship between the As quadrupole frequency and T c in 
oxy-arsenides.— 



IV. RELATIONSHIP WITH THE IRON 
SUPERCONDUCTORS 

To clearly show the relation of TlFe2Se2 with the Fe- 
based superconductors, we performed non-spin-polarized 
calculations of the electronic structure. The band struc- 
ture, DOS, and Fermi surface are shown in Figs. |U [5] 
and [51 respectively. The main features of band structure 
show close similarity with that of FeSe 18 and the Fe-based 
arsenidesi 36 ' 37 In particular, the states near Ep are domi- 
nated by Fe d states (mainly derived from xz, yz orbitals) 
forming a pseudogap at the electron count of 6 (rather 
than 4 in the tetrahedral crystal field scheme) , with some 
mixture with Se p states. Tl occurs as Tl + , as mentioned, 
with the Tl s states in the energy range ~ -8 eV to -4 
eV and the p derived valence states are above Ep. The 
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FIG. 3: (color online) Calculated total and projected elec- 
tronic DOS of TlFe2Se2 with the checkerboard antiferromag- 
netic order, using the LAPW method. For Fe, the minority 
half-filled (spin-up) and majority filled (spin-down) 3d states 
are shown separately. 
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FIG. 4: (color online) Calculated non-spin-polarized band 
structure of TlFe2Se2 using the relaxed zse = 0.342. 



resulting doping relative to FeSe (0.5 e/Fe) pushes Ep 
up from the steep lower edge of the pseudogap for FeSei 8 . 
to the bottom, leading to significantly decreased N(Ep) 
of 0.92 eV -1 per Fe. As a result, the hole Fermi sur- 
faces at the zone center in other Fe-based materials have 
completely disappeared and the electron sections at zone 
corner greatly expand, as shown in Fig. [SJ This destroys 
the strong nesting between two-dimensional electron and 
hole Fermi surfaces in more lightly doped Fe-based ma- 
terials and thus the SDW antiferromagnetic instability is 
suppressed. 

As in the other Fe-based superconducting materials^ 
the non-spin-polarized electronic structure of TlFe2Se2 
also shows great sensitivity to the Se height. For exam- 
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FIG. 5: (color online) Calculated non-spin-polarized total and 
projected electronic DOS of TlFe2Se2. The values are on a per 
formula unit basis. Note that projections are onto the LAPW 
spheres, thus the absolute contributions of the Tl-s and Se- 
p states are underestimated owing to their more extended 
orbitals which extend significantly into the interstitial region. 
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FIG. 6: (color online) Calculated non-spin-polarized Fermi 
surface of TlFe2Se2. The right panels are top views along 
Z — r direction. Here Fermi surfaces were mapped in the 
body-centered tetragonal Brillouin zone, and thus X point 
corresponds to the zone corner M point in the tetragonal 
Brillouin zone. 



pie, if the experimental z§ c = 0.357 was used for calcu- 
lations, N(Ep) increases to 2.0 eV _1 per Fe though E-p 
still lies in the bottom of pseudogap, and furthermore 
very small cylindrical hole Fermi surfaces appear at the 
zone center. This increase in N(Ep) with z§ c provides 
an explanation for the fact that, although not the ground 
state, an itinerant ferromagnetic state is lower in energy 
than the non-magnetic state even though the non-spin- 
polarized N(Ep) (with the non-spin-polarized zs e ) is well 
below the Stoner criterion for itinerant magnetism. 

The calculated V zz of Fe and Se are l.OxlO 21 V/m 2 
and -4.4xl0 21 V/m 2 with experimental z Sc and -0.4xl0 21 
V/m 2 , and -5.2xl0 21 V/m 2 with relaxed zs e , respectively. 
This illustrates the strong sensitivity of the EFG to the 
Se position, as mentioned above. Strong sensitivity of 



the EFG was also found for As in LaFeAsO by density 
functional calculations^ 



V. SUMMARY AND DISCUSSION 

To summarize, based on density functional calcula- 
tions, we study magnetic properties, electronic structure 
and relation to superconductivity for TlFe2Se2 occur- 
ring as the ThCr2Si2 structure. We find the nearest- 
neighbor checkerboard antifcrromagnetism is the ground 
state, consistent with the available experimental data. 
The non-spin-polarized electronic structure of TlFe2Se2 
shows close similarity with the Fe-based superconduc- 
tors. Relative to those materials, Tl + is an electron 
dopant, donating 0.5 additional carrier per Fe relative 
to Fe-Se layers. This over-doping significantly enlarges 
the electron sections of Fermi surface at zone corner and 
eliminates hole sections at zone center, and thereby com- 
pletely destroys the Fermi nesting and thus suppresses 
the SDW antiferromagnetic instability. We note that the 
compound crystallizes in the ThCr2Si2 structure, with a 
larger interlayer separation than FeSe. This might be fa- 
vorable to superconductivity by empirically considering 
that the highest T c so far is realized in the ZrCuSiAs- 
type compounds, i.e. the family with the largest inter- 
layer space^i 4 - In any case, TlFe2Se2 is not reported to 
have excess Fe in interstitial sites whose moments may in- 
teract with superconducting Fe-Se planes and cause pair 
breaking. However, some Fe vacancies might form as 
the intrinsic defects in this systemj 27 i 39 One natural way 
forward is via Tl deficiency, which would reduce the over- 
doping. 

In order to evaluate this possibility, we calculated the 
formation energy of a Tl vacancy, using a supercell as 
described above. The calculated energy of -0.14 eV, indi- 
cates metastability of the compound at low temperature, 
similar to PbZr0 3 rMI and that abundant Tl vacancies 
are very prone to form in this material. Thus, it may well 
be possible to produce Tl deficient material by modifying 
the growth conditions^ 2 - Owing to absence of excess Fe 
that cause pair breaking and suppress superconductivity 
in FeSe, if appropriate Tl deficiency introduced, this ma- 
terial (TLjFe2Se2) might be a good candidate for higher 
T c in chalcogenide family. This structure type may also 
provide an avenue for producing superconductivity in Fe- 
Te compounds, e.g. Tl a: Fe2Te2. 

Even in the over-doped regime represented by 
TlFe2Se2, where the Fermi level has dropped into the bot- 
tom of the pseudogap, resulting in relatively low N(Ep), 
we still found competing itinerant ferromagnetic and 
checkerboard antiferromagnetic states. In this regard, 
recent NMR experiments have shown evidence for pseu- 
dogap behavior in electron doped compounds (note this 
is the NMR pseudogap, not the pseudogap in the DOS), 
perhaps related to that in cuprates. 43 ' 44 However, the 
behavior is different. In cuprates the pseudogap appears 
to be closely associated with the magnetic ordering of 
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the undoped phases, and is seen most strongly in the 
underdoped regime. In contrast, the Fe-based supercon- 
ductors show the pseudogap in the electron over-doped 
regime and not in the underdoped regimei 45 ' 46 However, 
while Fe 2+ is a common valence for Fe, Fe 1+ is not and 
perhaps for this reason heavily over-doping beyond the 
superconducting regime while maintaining high sample 
quality has been difficult. Within an itinerant picture 
one may expect competition of SDW order and other 
magnetic (i.e. checkerboard) order in the Fe supercon- 
ductors, which gives way to checkerboard order when the 
hole Fermi surface, and therefore the nesting is destroyed 
by over-doping. Within such a scenario, and considering 
the experimental results, the pseudogap may be associ- 
ated with incipient checkerboard order. In this regard, 
the chemical stability of TlFe2Se2 may provide a very 
useful window into the over-doped regime, especially if it 
can be chemically doped with holes, e.g. by Tl deficiency, 
as discussed above. 

In fact, such such studies might be particularly illumi- 
nating considering the phase diagram that is suggested 
by the present results in conjunction with the known 
phase diagrams of the electron doped superconductors, 
first principles results for them and NMR data. In par- 
ticular, without doping the materials show a spin density 
wave, though with strongly reduced moments compared 
to standard density functional calculations. Such calcu- 
lations do not include renormalization due to spin fluc- 
tuations - an effect that is small in most magnetic ma- 
terials, such as Fe metal, but is apparently large in these 
materials. Electron (or hole) doping destroys the SDW, 
presumably by reducing the nesting between the hole and 
electron Fermi surfaces, in favor of a metallic state with 
evidence for strong spin fluctuations and superconductiv- 
ity possibly related to spin fluctuations connected with 
the SDW. In this regime, there also appears to be a strong 
renormalization of magnetism, based on comparison of 
density functional results and experiment j 34 i 47 In gen- 
eral such renormalized states occur near quantum criti- 
cal points or in situations where there is a competition 
between different magnetic orders. Considering the large 



composition range over which this behavior is seen, com- 
petition with another magnetic state may be crucial. The 
experimental observation that TlFe2Se2 is magnetic, and 
the likely checkerboard order suggests that the competing 
state in the superconducting phases might be the nearest 
neighbor antifcrromagnetic state. This state, although 
itinerant, is not driven by Fermi surface nesting, and 
therefore nearness to it would be expected to yield spin 
fluctuations over a broader region of the Brillouin zone 
than nearness to a nesting driven SDW. While this would 
make them more effective in renormalizing the SDW in- 
stability due to the larger phase space for competing fluc- 
tuations, it would also make them harder to observe by 
inelastic neutron scattering as compared to the nesting 
related fluctuations. Nonetheless, it would be of great 
interest first of all to confirm experimentally that the 
ground state of TlFe2Se2 is the checkerboard state, and 
secondly to study the overdoped regime, perhaps using 
this chemistry, to relate spin fluctuations from neutron 
scattering and the NMR pseudogap. Furthermore, it will 
be of considerable interest to study the interplay between 
superconductivity and the onset of checkerboard order, 
if this in fact is established at high doping. While the 
SDW and superconductivity mediated by nesting related 
spin fluctuations would compete for the same electrons 
at the Fermi surface, this is not necessarily the case for 
the checkerboard order, and so some difference may be 
anticipated from the behavior at the underdoped antifer- 
romagnetic (SDW) - superconducting boundary, perhaps 
including coexistence of the two orders over some com- 
position range. 
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